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Abstract 


This thesis attempts to develop a dynamic stability index of a general manipulator. It also 
deals with the analysis of dynamics and vibration of a 6-UPS Stewart platform. Finally, it 
presents a vibration control strategy for the same. 

The dynamic formulation for a 6-UPS Stewart platform follows the Newton-Euler ap- 
proach. Leg stiffness, force and torque due to viscous friction at the joints, inertia and gravity 
effects are considered in the model. 

The response of the platform, subjected to disturbances at different frequencies, has been 
studied. A PD control strategy has then been proposed for controlling both free and forced 
vibrations and is applied to isolate the top platform from base vibrations. 
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Chapter 1 


Introduction 

1.1 Introduction 

Technical applications resort to several instruments and machineries for achieving a concen- 
trated objective. The sensitivity of the paraphernalia is a critical point requiring careful 
attention. The key natural phenomenon against which these need safeguarding is vibrations. 
Objects which require or generate motion are invariably subject to detrimental effects of vi- 
bration where a huge analytical effort is spent for eliminating or reducing it. Let us take the 
all-too-familiar example of a naval tracking camera. These are used for reconnaissance and 
target locking in ships. They are subjected to ship’s motion viz. roll, pitch, etc. Similarly, 
reading of sensitive measuring instruments axe affected by the vibrations of the shop-floor, 
caused by other machines. So, it is important to isolate these devices from surrounding dis- 
turbances. For active vibration isolation, the instruments can be put on the end-effector of a 
manipulator. The manipulator may be serial, parallel or hybrid in nature. 

The advantage of serial manipulator is that its workspace is large. But due to cantilever 
nature, it is not very rigid and chatter is high. The links closer to base for a serial manipulator 
are subjected to the weight of the succeeding links in addition to the end-effector load. As 
a result, links are heavy for serial manipulators, rendering them unsuitable for high speed 
operations. High load and inertia also increase the requirement of power to drive them. 
Hence powerful motors are required for serial manipulators. 

In parallel manipulators, the end-effector is supported and articulated by independent 
closed-loop kinematic chains. As the load is divided between the supporting links, each 
link has less inertia and motor power requirement is less stringent. Closed structural form 
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of parallel manipulators reduces chatter and deflections. Because of this reason and also 
accentuated by the fact that the inertia is less, they are suitable for high speed operations. 
Besides, for parallel manipulators, the end-effector can be positioned more precisely because 
the deflections are quite small. On the other hand, small workspace associated with parallel 
structures forms another issue in this fleld. 

In some sense, the net motion at the end-effector, for a serial manipulator, is the summation 
of the individual link motions. On the other hand, end-effector force in parallel manipulators 
is the sum of individual link forces. Hybrid manipulators incorporate the advantages of both 
serial and parallel manipulators. They are made by combining serial structures in parallel or 
by parallel structures in serial manner or combinations thereof. 

1.2 Comparison Between Serial and Parallel Manipula- 
tors 

Each kinematic chain of parallel manipulators is constructed as a serial chain and hence they 
impose similar restrictions as serial manipulators. Therefore complex kinematic chain parallel 
manipulators cannot be compared in the following sense. 

The most significant differences between a serial and a fully parallel manipulator are the 
following: 

1. Degrees of freedom (DOF) of a serial manipulator is given by 

DOF = ^/, 

i=l 

where fi is the DOF of the z-th joint of the manipulator with a total number of j joints. 
For parallel manipulators, the expression of DOF is 

j 

DOF = m{n- j -1) ->rY^fi 

i=l 

where n is the number of links and j is the number of joints. The DOF of a rigid body 
in the ambient space is denoted by m. It is 3 for planar and 6 for spatial manipulators, 
respectively. 
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2. Direct kinematics involves finding the position and orientation of the end-effector for 
given joint angles and is very simple and unique for serial manipulators. They can 
easily be obtained by using Denavit-Hartenberg parameters (D-H parameters) and ho- 
mogeneous transformation matrices. For parallel manipulators direct kinematics is very 
complicated and requires a system of non-lineax coupled algebraic equations to be solved. 
Moreover, direct kinematics is not unique for parallel manipulators. 

Prom the concept of serial-parallel duality, inverse kinematics is complicated and have 
multiple solutions for serial manipulators. But it is simpler and often gives unique 
solution for parallel manipulators. 

3. The relation between joint velocities (©) and end-effector velocity (V) for serial manip- 
ulators is given by, 

V = J(©)© 

where J(0) is the Jacobian matrix. Its dual in parallel manipulator is force transfor- 
mation matrix [H(©)]. It relates joint torques (t) with the force supported by the 
manipulator (F), in following manners, 

r=H(©)F . 

At singular configurations, the Jacobian matrix of serial manipulators becomes singu- 
lar. Hence the end-effector loses one or more degrees-of-freedom making a part of the 
workspace inaccessible to the end-effector. On the contrary, at singular configurations, 
force transformation matrix of a parallel manipulator becomes singular. Hence it cannot 
support arbitrary load and gains one or more extra degrees of freedom. 

4. Redundancy is incorporated into a serial manipulator to make it more dexterous. This 
increases the workspace and also gives it the ability to avoid obstacles to reach out for 
a particular TSL (task-space location). In spite of this, redundancy cannot remove the 
singularity of serial manipulators. Singularity only gets redistributed. Detailed discus- 
sions on redundancy can be found in works of Nakamura [38]. In parallel manipulators, 
redundancy improves its force supporting capabilities at the cost of workspace volume. 
Unlike serial manipulators, singularity of parallel manipulators can be reduced by incor- 
porating redundancy. In certain cases, it can be successfully eliminated also. Dasgupta 
and Mrithyunjaya [6] have discussed it in detail. 
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5. Dynamics of serial manipulators are expressed in terms ordinary differential equations. 
In case of parallel manipulator these relations are exacerbated in complexity by con- 
straints as they quite often turn out to be differential algebraic equations. 

6. Control of serial manipulator is relatively simpler and can be done by proportional plus 
derivative (PD) control scheme. Control of parallel manipulator is much more complex 
due to nonlinear coupling. 

1.3 Stewart Platform 

Stewart platforms are most widely studied among the parallel manipulators. In general, Stew- 
art platform can be described as two rigid bodies, namely base and top platforms, connected 
by six linearly actuated legs. Having six parallel-actuated simple kinematic chains, its DOF is 
also six. Connections between legs and platforms may be both spherical joints (called 6-SPS 
Stewart platform (Fig (1.1)) or one spherical and other universal joint (called 6-UPS Stewart 
platform (Fig (1.2)). 
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Figure 1.2: The 6 - UPS Stewart Platform 

1.4 Literature Review 

The Stewart platform derives its name from that of Stewart [35], who designed a mechanism 
for generating general motion for flight simulation. Gough suggested the use of six linear 
actuators in- parallel to actuate the platform [16]. McCallion and Truong first designed a 
Stewart platform for mechanical workstation [25] . They also performed some theoretical and 
numerical studies on mobility analysis and iterative solution of direct kinematics. Detailed 
study on workspace section and kinematics of Stewart platform were done by Merlet [29], 
Fichter [9], Yang and Lee [37]. 

Closed form dynamic equations of Stewart platform assuming rigid legs and stationary 
base were derived following three approaches viz. Newton-Euler method [5, 4], Lagrange 
method [19, 11] and virtual work principle [40]. 

Lee et al. [21] used Stewart platform for control purposes assuming legs to be rigid and 
base as fixed. Ren et al. [13] and Cheng et al. [2] used Stewart platform for isolating wind 
excited vibrations on a telescope mounted on a hanging cable car. They accounted for the 
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movement of the base but assumed components of the Stewait platform to be rigid. With the 
same assumptions, Kim and Lee [17] developed the mathematical model of Stewart platform 
for isolatin.g a slender bar from the base excitation. 

Shuguang and Schimmels [33] studied the behaviour of simple springs connected in parallel 
to a rigid body. Selig and Ding [32] developed a mathematical model of Stewart platform 
for the study of vibration, with simplified assumptions i.e. massless legs, frictionless passive 
joints. Lee and Geng [20] derived the dynamic model of Stewart platform by considering 
flexibility of its legs and assumed the base to be static. For micro-scale operations, Mclnroy 
[26] used flexure jointed hexapods. There he did not consider the inertia of the legs and 
assumed the base to be static. When we are interested in a particular position and it is 
pre-decided, the Stewart platform can be designed for this particular position. Geng and 
Haynes [12] have deflned a special architecture of Stewart platform, which they termed as 
‘cubic configuration’, where there is no coupling between the DOF of the Stewart platform. 
Hence, vibration control turns out to be simple at this configuration. Once it is fabricated, it 
can be used for this particular position and orientation only. 

Basic study of singularities in parallel manipulators was done by Gosselin and Angeles 
[15]. They classified the singularities in 3 types. ‘Singularity of 1st kind’ (motion singular- 
ity), ‘singularity of 2nd kind’ (force singularity), ‘singularity of 3rd kind’ (both singularities 
simultaneously). Stewart platform can have only ‘singularity of 2nd kind’. Later Ma and 
Angeles [23] showed that a highly symmetric Stewart platform (e.g. Stewart platform with 
regular hexagonal base and top platforms) is always singular. He called it ‘architecture singu- 
larity’. Gosselin [14] also studied singularity associated with high condition number of force 
transformation matrix. Bhattacharya et al. [1] and Dasgupta et al. [7] proposed strategies 
for singularity free path planning. They used a high condition number of force transforma- 
tion matrix as the representation of singularity. In all these cases, singularity studies have 
taken into account the architecture and configuration (i.e. kinematic parameters) of Stewart 
platform. The dynamic parameters (i.e. mass, inertia, stiffness etc.) were not taken into 
account. 

Yoshikawa [39] proposed ‘dynamic manipulability ellipsoid’ (DME) to take care of the 
dynamic parameters of manipulators. It signifies, at a particular configmration how arbitrary 
change in acceleration of the end-effector can be obtained with joint driving forces. Though it 
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takes care of the inertia parameters of manipulator, it does not consider the stiffness. Ma and 
Angeles [22] have introduced ‘dynamic conditioning index’ (DCI) and used it for designing 
manipulator for trajectory planning [24]. A low value of DCI. signifies that the inertia matrix 
is close to ideal isotropic generalized inertia matrix. For an ideal isotropic generalized inertia 
matrix, inertia torques are completely decoupled, wherein they are easy to control. 

1.5 Motivation and Scope of the Thesis 

The previous section indicates that though kinematic and architectural singularities have been 
studied in great detail, effects of dynamic parameters have not received adequate attention. 
If the manipulator is used for control, it has to be dynamically stable. Hence there should be 
some well defined index for dynamic stability. Using Stewart platform as a representative, a 
dynamic stability index has been proposed in the 1st part of the present work. 

The literature survey also indicates that, though detailed studies of dynamics of Stewart 
platform have been conducted considering fixed base and rigid legs, effect of moving base and 
flexible legs were not studied in great detail. In the second part of this work, completely 
general closed form dynamic equation of Stewart platform, following Newton-Euler approach, 
have been developed, considering the base movement and also accounting for the stiffness of 
the legs. 

Subsequent to the above steps, we have employed PD control law to isolate the top platform 
from vibrations being transmitted from the base. Simulations have been performed using the 
closed form dynamic equations. 

1.6 Thesis Organizations 

In the next chapter, we develop a stability index which will take care of both kinematics and 
dynamics of manipulators. Comparisons of this with general singularity measure has been 
made through examples. 

In chapter 3, completely general dynamic equations of Stewart platform have been devel- 
oped using the Newton-Euler approach. In chapter 4, the stability index, developed in chapter 
2, has been validated through numerical simulations. 

In chapter 5, a linearized model, for gain selection, has been developed. For chosen values 
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of gains a PD control law is implemented to isolate the top platform from base vibrations and 
simulation results are included. 

A summary of the work and future scopes have been discussed in chapter 6. 
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Chapter 2 

Dynamic Stability Index 

2.1 Introduction 

FVom literature survey in Chapter 1, we observe that the effects of the dynamic parameters 
of a manipulator (i.e. mass, inertia, stiffness etc.) have not been taken into account for 
describing its stability. The natural frequencies, determined from linearized dynamic model 
of a manipulator, incorporate the contributions of all kinematic and dynamic parameters of 
the manipulator. Therefore, the lowest natural frequency may very well represent the stability 
of a manipulator. A low value of it signifies that the manipulator is tending towards neutral 
stability. Hence, it easily gets perturbed from desired position under small disturbances. These 
have been explained in this chapter by taking 6-UPS Stewart platform as an example. 

In the next section, first we linearize the dynamic model to find out the natural frequencies. 
Then the sensitivity of the lowest natural frequency to dynamic parameters has been studied, 
taking numerical examples. Conclusions have been drawn from the simulation results. 

2.2 Linearization of Dynamic Model of Stewart Plat- 
form for Free Vibration 

We have assumed, 

1. The base and top platforms are rigid. 

2. Legs have lumped constant stiffness (due to hydraulic actuators, material compliance 
etc.) in axial direction. But they have no rotational stiffness. 
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The closed form dynamic equation of 6-UPS Stewart platform with fixed base in joint 
space [5] is 

, ( 2 . 1 ) 

where 

J = combined inertia of leg and moving platform. 

H = input-output force transformation matrix. 

rj = force due to gravity, Coriolis acceleration and viscous friction at joints. 

P = force, generated by actuators at the prismatic joint of legs or due to stiffness of legs, 

u = centripetal acceleration of joining point of the moving platform and leg. 

5ft = rotation matrix specifying the orientation of the moving platform [3]. 

Fext = resultant of external forces acting on the moving platform. 

Meii = resultant of external moments acting on the moving platform. 

We are now interested only in the platform response to initial disturbances, so we can take 
external loads to be zero, without losing any generality. So, we get Fext = Mext = 0. For free 
vibration, magnitude of the force in i-th leg is 

Fi = {KieMLi - Loi) = {Kteg)i5Li , (2.2) 

where 

Loi = equihbrium length of i-th leg, 

Li = current length of i-th leg. 

For free vibration, rj and u are negligible compared to other terms of Eqn. (2.1). So 
combining Eqn. (2.1) and Eqn. (2.2) we get 

= -K6L , (2.3) 

where 

K = dmg{{Kieg)i). 
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Here J and H both depend on the constant properties of the platforms and legs (e.g. 
position of platform and base joints, mass and centre of mass of moving platform and legs 
etc.) as well as position and orientation (t and 6) of moving platform. So, natural frequencies 
and mode shapes of Stewart platform also depend on position and orientation. For a particular 
position and orientation of a given Stewart platform , both J and H are constants. Hence 
is constant for that position and orientation. We denote as Jo. So 

Eqn. (2.3) becomes, 

JoL = -K5L. (2.4) 

Assuming to be a solution of the Eqn. (2.4), we get the a generalized eigenvalue problem 
[31] 

u^JoL = KL. (2.5) 

Solving Eqn. (2.5) we get six natural frequencies and six mode shapes [28, 27, 36] for a 
particular position and orientation of the Stewart platform. 

The lowest natural frequency, calculated from Eqn. (2.5), is an index for dynamic stability 
of the manipulator. When the value of the lowest natural frequency is zero, the manipulator 
becomes neutrally stable, and hence, it cannot restore its position and orientation if it is per- 
turbed by any amoimt. A very low value also signifies dynamic instability of the manipulator. 
Therefore, the lowest natural frequency plays the same role in dynamic stability as does the 
condition number of the force transformation matrix in statics. High value of the condition 
number of the force transformation matrix signifies the static singularity of the manipulator. 

2.3 Simulation 

The generalized eigenvalue problem (Eqn. 2.5) has been solved in MATLAB, taking two 
sets of numerical values for two different designs of Stewart platform. Between them, one is 
completely general [5] and other one is designed for minimum condition number of H [8]. The 
variation of the lowest natural frequency with condition number of H for different mass of 
the moving platform (M) is shown in Figs. (2.1) - (2.4). Suitable stiffnesses of the legs for 
both the cases and dynamic parameters for second case have been assumed. The values of the 
parameters are given in the appendix A. 

For different position and orientation of Stewart platform, though the condition number 
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of H may be same, the lowest natural frequency may be different. This can be seen from the 
plots also. 

1. Example: 1 

(a) Case: 1 M = 40 kg, (b) Case: 2 M = 4000 kg. 

2. Example: 2 

(a) Case: 1 M = 4 kg. (b) Case: 2 M = 400 kg. 

In Figs. (2.1) - (2.4) we can see that the lowest natural frequency of Stewart platform 
decreases very sharply with increase of condition number of H. As H tends towards ill- 
conditioning, the platform cannot support load coming from certain directions. Therefore, it 
gains one or more degrees of freedom. Very low value of the lowest natural frequency also 
signifies the instability of the moving platform. 

But, at the same condition number of H, at different poses (different position and orien- 
tation), we see a lot of variation of the lowest natural frequency. Here, condition number of 
H alone cannot point out which position and orientation is more stable. 

With the same kinematic parameters, but with different mass of the moving platform, 
variation of the lowest natural frequency is significant. With the change of mass of the top 
platform, drastic change of the lowest natural frequency of the test manipulator J, is clear from 
Fig. (2.1) and Fig. (2.2). Similarly, the variation of the lowest natural frequency of the test 
manipulator //, can be seen in Fig. (2.3) and Fig. (2.4). In both the cases, the lowest natural 
frequency decreases sharply with increasing mass of the platform, which signifies increasing 
dynamic instability of the moving platform. This effect also cannot be explained with the 
condition number of H. 

2.4 Conclusions 

Prom the plots and the discussion in the previous section we can make following conclusions. 

1 . Where dynamic stability is more important, design criterion can be maximization of the 
lowest natural frequency. When we are interested in different positions and orientations. 
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instead of a particular one, i.e. if the manipulator is going to be used in a path planning 
and tracking purpose, the design criterion can be that the lowest natural frequency is 
higher than a lower bound over the entire path or it should be maximized over the path 
in a cumulative sense with some appropriate weights. Here, more stringent first criterion 
needs more computation resources and also requires robust optimization routine to solve 
them. Little relaxation of the requirements can be done by following the second criterion. 
But in this case, appropriate weights should be decided based on the criticality of the 
operation. Over the whole path if there axe some poses where manipulator needs to be 
made more stable from the appMcation point of view, weights should be high for these 
poses. 

2. The natural frequencies are also important for choosing the gains of PD controller. Gains 
should be chosen in such a way so that the close loop poles of the corresponding linearized 
system (linearized model has been given in chapter 5) should not fall near the natural 
frequencies. If they fall near the natmral frequencies, the controlled system will vibrate 
more than the uncontrolled one. 

3. Prom the natural frequencies, we get a rough estimation of the samphng time (At) for 

controller. For proper control, sampling time should be At = — [10]. Where, 

is the highest natural frequency. 
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Chapter 3 

Closed Form Dynamic Equation of 
Stewart Platform with Base Motions 
and Flexible Legs 

3.1 Introduction 

Dynamics of Stewart platform is quite complicated due to its closed chain construction and six 
degrees of freedom. Prom literature survey, we find that a good amount of research has been 
pursued to find out closed form dynamic equations of Stewart platform, following both the 
Newton-Euler and the Lagrangian approaches. But the effects of base motion and flexibility 
of the legs have not been throughly explored. Use of a completely general dyneimic model will 
give better performance of the controller. In this chapter, a generalized dynamic model of 
a 6-UPS Stewart platform following the Newton-Euler approach has been developed in the 
same way as by Dasgupta and Mrithyunjaya [4, 5]. Base motion, friction at the joints, inertia 
of legs and platform, flexibility of the legs are also considered in the model. Leg stiffness, 
present only in axial direction, is assumed to be lumped. 

3.2 Closed Form Dynamic Equations for 6-UPS Stewart 
Platform 

To develop closed form dynamic equations, first kinematics and dynamics of the legs have 
been studied. Then the force, applied on the top platform by each leg, and the external loads 
have been combined to formulate the dynamic equations of the top platform in task space. 
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While deriving the equations of motion for a leg, the suffix i, denoting a general leg, is omitted 
to reduce the cluttering of the equations with too many indices. Prom the context, it is quite 
clear that this equation is for a general leg. At the time of writing equations for. the platform, 
the suffix i has been incorporated. 


X 


Figure 3.1: Positions of Reference Frames 




The detailed diagram depicting different reference frames and a schematic leg is shown in 
Fig. (3.1). The symbols used to develop the model are explained below: 

G = global inertial reference frame. 

B = base frame (attached to the reference point of the base platform) . 

P = platform frame (attached to the reference point of the top platform). 

D, U = Frames, attached to lower and upper part of the legs, respectively. 
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tb = position of the frame B with respect to frame G. 

tp = position of the frame P with respect to frame G. 

3^6 = rotation matrix (orientation of the frame B w.r.t. frame G). 

= rotation matrix (orientation of the frame P w.r.t. frame G). 

bj = i-th base point in the base frame. 

Pi = i-th platform point in the platform frame. 

<Jbi — S^ibj. 
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5^Pi. 

= linear velocity of the base platform w.r.t. frame G. 

= linear velocity of the top platform w.r.t. frame G. 

= angular velocity of the base platform w.r.t. frame G. 

= angular velocity of the top platform w.r.t. frame G. 

= linear acceleration of the base platform w.r.t. frame G. 

= linear acceleration of the top platform w.r.t. frame G. 

= angular acceleration of the base platform w.r.t. frame G. 

= angular acceleration of the top platform w.r.t. frame G. 

= centre of gravity of the top platform in platform frame. 

= mass of the top platform including payload. 

= moment of inertia of the top platform including payload in platform frame. 
= external force on the top platform in the platform frame. 

= external moment on the top platform in the platform frame. 

= stationary axis of the universal joint at i-th leg. 

= rotation matrix (orientation) of i-th leg w.r.t. frame G. 

= masses of the lower and upper part of i-th leg. 

= centre of gravity of the lower and upper part of i-th leg in local frames. 

= Moments of inertia of lower and upper part of i-th leg in local frames. 

= acceleration due to gravity. 
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Cu, Cp, Cs = coefficients of friction in the universal, prismatic and spherical 
joints respectively. 

Fi = magnitude of input force at i-th leg. 

Sj = vector along i-th leg. 

3.2.1 Velocity Analysis of i-th. Leg 

Vector along i-th leg: 

S = (qp + tp) - (qb + 1<,). 

Hence leg length: 

L=|1S1|. 

Unit vector along the leg: 

s = S/L. 

T=lx y z], 

where 

k = and x = s; y = (k x s)/ || k x s ||; z=xxy. 

Centres of gravity of lower and upper part w.r.t. global reference frame: 

Td = Trdo, 

r„ = T(v-|-r,»o), 

where 

v = [L 0 Of. 

Moments of inertia of the lower and upper part w.r.t. global reference frame: 

Id = TldoT^, 

= T[I„o + m„L2diag(0,l,l)]T^. 

Relative velocity between the base point and platform point: 

S = (ojp X qp -f tp) - {(jJb X q;, tb). 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


(3.5) 

(3.6) 


(3.7) 

(3.8) 


(3.9) 
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Hence relative sliding velocity between upper and lower part of the leg: 


L = s-S. 


(3.10) 


Angular velocity of the leg: 


W = sxS/L. 


(3.11) 


3.2.2 Acceleration Analysis of t-th Leg 

The relative acceleration between the base point and platform point: 

S = [tp + cvp X qp + Wp X {u)p X qp)] - [t{, + at x q^ + 0^6 x (w* x qt)]. (3.12) 


Let us group the velocity dependent terms in Ui and the acceleration dependent terms in Up 
and at, as 


Ui = cOpX {uip X qp) - wt X {iOb X qt), (3.13) 

B-p = tp + Op X qp, and (3.14) 

at = t't + atxqt. (3.15) 


Hence, S becomes 


S — Up at "f” XJ 1 . 


(3.16) 


S can be expressed in terms of leg quantities like relative sliding acceleration, angular velocity 
and angular acceleration as, 


S = Ls + W X (W xS) +2W X Ls + A X S. (3.17) 


Presence of a universal joint, at one end of the leg, inhibits its rotation about its own axis 
along the length. Therefore, W • s and A • s are zero. Equating S • s from equations (3.16) 
and (3.17) and using the relation W • s = 0, we get 

L = s • (ap - at) -f u (3.18) 


where 

u = s- lJi + y(S-Ls)\ 
Lj 


(3.19) 
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Similarly equating s x S from both equations (3.16) and (3.17) along with the relations 
W -8 = 0 and A • s = 0 we get 

A = isx(ap-ab) + U2 "(3.20) 

where 

U 2 = |(s X Ui -2LW). (3.21) 

L/ 

Accelerations of centre of gravity of lower and upper parts of the leg are respectively, 

ad = [tfe + Q!6 X Qb + Wj, X (cdb X qb)] + A X Fd + W x (W x rj), (3.22) 

au = [t'b + Ckb X qb + Wb X (wb x qb)] + Ls + A x 

+W X (W X r„) + 2LW x s. (3.23) 

Simplifying and using Eqn. (3.20), we get 

ad = ^(s X ap) X Td + [ab - ^(s x ab) x id] + U 3 (3.24) 

= (s • ap)s + y(s x ap) X r„ + [ab - y(s X ab) X Fu - (s •ab)s] +U 4 (3.25) 

Jb L/ 

where 

U3 = U2 X Fd + W X (W X Fd) + Wb X (cdb X qb) (3.26) 

U4 = ws + U2 X Fu 4- W X (W X F„) + 2LW x s + Wb x (wb x qb) (3.27) 

3.2.3 Dynamics of t-th Leg 

If Fs is the force acting on the leg at the spherical joint, from moment equilibrium of the 
entire leg, we get 

-rridTd X ad - rriuTu x a„ + {mdTd + muVu) x g - (Id + L) A 

-W X (Id + I„)W + M„s + S X F, - C„W - f = 0 (3.28) 


or, 

M„s + SxF, = C , (3.29) 

where 

C = rUdVd X ad + m„r„ x - (mdrd + rriutu) x g 

+(Id + I„)A + Wx(Id + I„)W+C„W + f. _ (3.30) 
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The expression of f, the moment acting at the spherical joint due to viscous friction, is 

f=C,(W-a;p). (3.31) 

Taking cross product of both the sides of Eqn. (3.29) with s and rearranging, we get 

F, = (s • F,)s + |c X s. (3.32) 

Lj 

Prom the force-equilibrium along the axis of the leg, 

F + s • F^ - CpL + rUuS • g = m^s • a^. (3.33) 

Substituting s • Fj from Eqn. (3.33) into Eqn. (3.32) and using Eqns. (3.20), (3.26), (3.27) 
and (3.30) we get, 

Fs = mu[s • ap + i{(s • ru)s • ap- r„- a^} - i{(s • r„)s • afc - Fu • ab}]s 

-js X [mdTd X {y(s X Up) X r<i + (ab - y(s x ab) x r<i) 

Ju Lj 

■^TriyTu X {(s • ap)s + ^(s x ap) x Tu + (ab - ^(s x Ub) x - (s • ab)s)} 
+y(Id-l-IOs X (ap-ab)]+V-sF (3.34) 

Lj 

where 

V = (muS • U 4 -I- Cpi - m^s • g)s — ys X U 5 (3.35) 

Lj 

and 

U5 = rridYd X U3 -b rriuTJu x U4 4- (Id + Iti)U2 

-l-W X (Id -b I„)W - {mdTd + rn„r„) x g -b C^W -b f. (3.36) 

Eqn. (3.34) can be written in compact form: 

F, = Qpap - Qbab + V - sF (3.37) 


where 

Qp 


[m„(l-b— ^) ^2 ps + 1 E3 


m. 


1 


L2 


^(srl + r.s^) 


^[md(s X rd)(s X Td)^ -bm„(s x r„)(s x r„)^ -b s(Id -b I„)s] 


(3.38) 


Q6 = 


s • r,. 


f rridrl + m^rl , rudr^ + m^r^ rnu,_T , ^ 

[”^u— ^-^2 J®® 1 ^3 - + s r„E3j 

-b^(rdS^ - s^rdEs) - ^[md(s x rd)(s x rd)’" 

-bmu(s X ru)(s X r^)^ + s(ld-bl„)s] 


(3.39) 
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Es is the 3x3 identity matrix and 

0 -S^ Sy 

s = 0 -Si 

Now substituting expressions for ap and from Eqns. (3.14) and (3.15), respectively, into 
Eqn. (3.37) and simplifying, we get the expression of force acting at spherical joint of z-th leg 
as 

(Fs)i = (Qpi'ip - Qbi'ib) - (QpiqpiCtp - QbiqbiQ;6) + Vi - SiFi. (3.40) 

3.2.4 Kinematics of Top platform 

Position vector of center of gravity of top platform in global reference frame is 

Rp = 3?pRop. (3.41) 

The acceleration of top platform in global reference frame is 

a = tp + Up X (u>p X Rp) + Op X Rp. (3-42) 

The moment of inertia of top platform including payload in global reference frame is 

I = 3fipIp5Rj. (3.43) 

3.2.5 Dynamics of Platform in Task-Space 

Prom force balance of top platform, we get 

6 

-Ma + Mg + SRpFext - ^(F.)i = 0. (3.44) 

2=1 

Substituting expressions of a and (F^)^ from Eqns. (3.42) and (3.40) into Eqn. (3.44) and 
simplifying we get 

6 ^ ^ ® 

(AfEs + ^ ^ Qp Jtp — (MRp + y ^ 

2=1 i=l i=l 

6 6 6 

+ + M{ujp X (u^p X Rp) ~ g} H- Vj = SiFi + ^pFext- (3.45) 

2=1 i=l i=l 
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Using Euler’s equation for platform about the platform reference point we get, 


— MRp X [tp + ttp X Rp + X (a;p X Rp)] + MRp x g — lap — UpX lup 

6 6 

“b ^ fi H” ^ CJp^ X (Qb^tp Qb^Qb^^^p) 

i=l i=l 

6 

— ^ ^ qp^ X {Qp^tp — Q^pi^piO'p + Vi — SiEi} = 0 . ( 3 . 46 ) 

i=l 

Again, substituting expressions of a and (F^)^ from Eqns. (3.42) and (3.40) into Eqn. (3.46) 
and simplifying, we get 


(AfRp + ^ ^ QpiQpi)^p “ QpjQbjtp 
i=l i=l 

6 

+[I + M(RlE, - R,Rj) - 

i=l 

6 

”t" ^ Q.PiQbiQ,biQ:p "I* X I^p 
i=l 

6 6 

H-MRp X [(a;p ■ Rp)a;p — g] + x Vi — f;) = ^ ^i)^i + S^pMeit- (3.47) 

i=l i=l 


Combining Eqns. (3.45) and (3.47) we get the closed form dynamic equations as 


Jp 



tb 

Oib 


T} — JJF -|- 


^pFext 

^pMext J ’ 


(3.48) 


where 




6 


'plat 


+EJ. 


^plat 


Jpi 


MEz -MRp 

MRp l + MiRlEz-RpRl) 

Qpi 

Qpi^ipi ~%>iQpi%i 


Ji, 


6 


EJb., 

i=l 


J6. = 


Qbi “Qbiqbi 

QpiQbi QpiQbiqbi 


J 
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6 

Vplat + ^ ^ 

1=1 

M{ujp X {ujp X Rp) - g} 

_ ojp X lujp + MRp X {(Wp • Rp)6jp - g} J ’ 

Vi 1 

. qpi X Vi - fi J ’ 

. qpi X Si qp2 X S2 qp3 X S3 qp^ x S4 qp^ x S5 qpg x Se J ’ 

[Fi F2 F3 F4 F5 Fe]^. 

The forces, coining from the base platform and exerted by the actuators, will be transmitted 
to the top platform through the legs only. As the legs are flexible, there is a small relative 
acceleration between any two points in the legs. Hence an inertia force also comes due to 
it. But for a Stewart platform, constructed for practical application (for vibration isolation, 
tracking, sensor application, etc.), the stiffness of the legs are enough to make this inertia 
force negligibly small. Therefore, the net force transmitted to the top platform through a leg 
is equal to the force generated due to the deformation of the legs. So, for z-th leg, 

Fi = {Kieg)i(Loi — Li) 


V = 

Vplat ~ 

Vi = 

H = 

F = 


where 


Loi = equilibrium length of z-th leg. 

Li = current length of z-th leg. 

3.3 Conclusions 

Using the Newton-Euler approach, the dynamic formulation of a 6-UPS Stewart platform 
has been carried out in the previous section. This model can be used for design and control 
purposes. For designing a Stewart platform some forces and moments, the expressions of 
which have not been given so far, are also necessary to be calculated to find the strength and 
deflection of the structure. Those can be easily found out from the forces and accelerations, 
already determined. 

Here we have assumed lumped stiflfness of the legs and neglected the inertia forces, coming 
on the legs due to their flexibility. Using finite element and finite difference method these 
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restrictions can be avoided. But the finite element models, being too computation intensive, 
would make simulations very slow. This hindrance also can be circumvented by exploiting 
the inherent parallel nature of the Stewart platform . The calculations for six legs can be 
done independent of each other, only with the information of base and top platform motions. 
Owing to this fact, use of parallel programing will drastically reduce the computational time. 
These can be explored in a future work in this area. 
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Chapter 4 

Validation of Linearized Model 


4.1 Introduction 

In chapter 2, we have used the lowest natural frequency as a measure of stability of a manip- 
ulator. Calculation of natural frequencies from the linearized dynamic model was explained 
by taking 6-UPS Stewart platform as an example. In this chapter with the help of numerical 
simulations we have justified that linearization. 

4.2 Numerical Simulations 

The dynamics of Stewart platform is expressed by a system of stiff, coupled, nonlinear, second 
order, ordinary differential equations [34, 18]. Till now, no analytical solution of these kind of 
equations has been found out. Taking position, orientation, linear and angular velocities as 
state variables, we express Eqn. (3.48) in state-space as 




z(7: 12) 



z = 

J-1 ^HF-t- 

3ftpFext 

— ’7 -h J6 


) 


L V 

_ 5ipMeit 


. 

) J 


Here z is the state vector, and is given by 

T 

^ ~ [ ^Vx ^Py ^Pz ^Px ^Py ^Pz ^Pz ^Py ^Pz ^Px ^Py ^Pz ] ’ ( 4 - 2 ) 

Numerical solution of Eqn. (4.1) have been performed in MATLAB for the test manipulator 
II, given in appendix A, for the top platform’s reference position, (0, 0, 0.085) and orientation, 
(0, 0, 0), subjected to base excitations at different frequencies. We get two distinct nature of 
responses of the top platform, depending on whether the excitation frequencies are equal to 
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Translation in Z direction Rotation about Z axis 


Figure 4.1: Displacement from Equilibrium Position Vs Time Plot at Base Excitation Frequency 
Equal to 2nd Natural Erequency, u = 749.715 Hz 
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Translation in Z direction Rotation about Z axis 


Figure 4.2: Displacement from Equilibrium Position Vs Time Plot at Base Excitation Frequency, 
uj = 550.136, Away from Natural Frequencies 
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the natural frequencies or away from it. Two typical plots one at the 2nd natural frequency 
(Fig. 4.1) and other at away from the natural frequencies (Fig. 4.2) have been shown. 

For different base excitation frequencies, we get the results shown in Fig. (4.1) and Fig. (4.2), 
depending upon whether the frequency matches with the natural frequencies or not. For dif- 
ferent base excitations if we plot the ratio of the response amplitude to the base excitation 
amplitude against the base excitation frequency, we get a plot like Fig. (4.3). Prom the dy- 
namic model, it is clear that the Stewart platform has six natural frequencies. (With flexible 
legs the platform actually it has infinite number of natural frequencies, due to approximations 
it comes to six.) Computations have been carried out for base excitation frequency up to a 
little above the second natural frequency. 

4.3 Conclusions 

It is clear from Figs: (4.1) and (4.3) that the Stewart platform exhibits resonance at the 
natural frequencies, calculated from the linearized dynamic model. Hence these frequencies 
can be considered as a representative of the system dynamics. If the lowest natural frequency 
is very small the manipulator tends to be unstable and fails to support loads or to track a 
trajectory. 
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Chapter 5 

Vibration Control 


5.1 Introduction 

Originally, the Stewart platform was developed for flight simulations. But, soon, it found its 
use in several applications, due to its six-DOF motion capabilities. High load carrying capac- 
ities and precise positioning capabilities give it a wide range of applicability from supporting 
missile launch pad to micro-scale sensor applications. One of the application is as a vibration 
isolator. Prom literature review we see, though many types of controllers have been designed 
assuming rigid body model of the Stewart platform, its flexibility is not taJken into account. 

In this chapter, a proportional plus derivative (PD) control strategy for isolating any 
device, mounted on a Stewart platform, from surrounding vibrations has been elaborated, 
considering the flexibility of the legs. For calculation of gains, a linearized model of the 
Stewart platform is developed. Then, using the PD controller, numerical simulations of free 
and forced vibration control have been performed to evaluate the system performance. 

5.2 Vibration Control Strategy 

We have used PD control in task space. We can get the displacements and velocities of the 
reference point of the top platform in task space in following ways. 

• We can directly measure the displacement of the platform reference point by camera 
sensor and get its velocity by using numerical differentiation [31, 18]. In this case, 
central difference scheme will be preferable. 

• We can obtain the acceleration of the platform reference point by accelerometer and get 
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displacement and velocity by numerical integration [31]. 

• We can measure the change of leg-length and its rate by means of instrumentation of the 
active joints. From that, we can solve the forward kinematics to find out displacement 
and velocity of the platform reference point in task space. (In general forward kinematics 
of parallel manipulators requires solution of a set of coupled nonlinear algebraic equations 
which is computation intensive. But, as the next point is always very close to the 
previous point so Newton-Raphson or other local methods [31] wiU very quickly converge 
if the previous point is taken as an initial guess for the numerical routines. Only for the 
initial point it will take considerable time.) 

• Any combination of the above methods. 



I CONTROLLER 


Figure 5.1: Control Diagram 


If there were six actuators which could apply Cartesian forces directly on the top platform 
to stabilize its motion, let u be the force vector apphed by them. Here u is a 6 x 1 vector 
and each component of it signifies the force given by each of those virtual actuators. Then 
the dynamic Eqn. (3.48) becomes 



+ 7] = HF + 




(5.1) 
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Let the force vector u be generated by a virtual PD controller having gains Kg. If 5z is the 
perturbation of the state vector from its equilibrium position Zq, then the force generated by 
the actuator is 


u = -K.gSz (5.2) 

where 


(5z = z(t) - z(0) = z[t) - zo 


The gain matrix is partitioned into two parts: proportional gains, Kp and velocity gains, Kt,. 
Therefore 

Kg = [Kp K„]. (5.3) 


Eqn. (5.1) can be written in state space as 



z(7: 12) 




r 0 1 

1 

+ 

_ SRpMeit _ 

-V + ^b 

1 1 

o- 


+ 



(5.4) 


Here, Jf, 



is not known a priori. But depending on the applications, a rough estimation 


of the maximum magnitudes of t(, and can be made. At the desired control position, we 
know the value of Jb. Hence a preliminary approximation of the whole term can be made, 
and denoted as d. 

The first term of the right hand side of Eqn. (5.4) can be expanded in Taylor’s series 
around equilibrium point (zq) as following manner 


First term = /(zq) + Adz -b Higher order terms (5.5) 

where 




z(7 : 12) 


A = Jacobian of 

+ 

1 

1 

3?pFeit 

_ JftpMext 

-r? + d^ 


The expression of Jacobian is developed in the appendix B. 

As zo is the equilibrium point, 

/(zo) = 0. (5.6) 
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Using Eqn. (5.2), the second term of Eqn. (5.4) can be written as 

K/z = -BKp(5z (5.7) 

where 



Prom Eqn. (5.5), Eqn. (5.6) and Eqn. (5.7), we can write Eqn. (5.4) up to first order approx- 
imation as 

z{t) = {A-BKg)5z{t). (5.8) 

Hence for no steady state error the eigen values of (A — BK^) have to be on the left half 
of s-plane. Now by any standard pole placement method we can get K^. (Refer Ogata [30], 
for a detailed discussion on pole placement methods.) 

The gains evaluated above would have stabilized the platform if the controlling forces 
instead of through legs were applied directly on the top platform by actuators. Hence if the 
force transformation matrix at the desired position is Hq, the actual gains are, 

K, = Ho-‘K,. (5.9) 

Now Kg consist of actual proportional and derivative gains 

Kg = [Kp K,]. (5.10) 

The control force generated by the actuators at the legs are 

u = -Kg5z. (5.11) 

As this control force, u is applied through leg actuators, they will also compress the legs and 
will be a part of leg force F along with the spring force generated due to base motions. 

5.3 Simulations 

Using the vibration control strategy, described in the previous section, numerical simulations 
of control of actual plant against free vibrations and vibration due to base motions have been 
made. Some of the results are presented in the next two sections. 
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5.3.1 Free Vibration Control 

For the simulation of free vibration, a small instantaneous displacement of the top platform 
has been given in absence of any other external disturbances. Numerical data of the platforms, 
given in the appendix A, have been used for simulations. 

All numerical values are in SI units. As the positions and velocities of the Stewart 
platform are strongly coupled, hence, general gain matrices have been used for controlling 
the vibration instead of diagonal gain matrices. These coupled gain matrices improve the 
controller performance. 

1. Case 1: Test manipulator I 


Desired position: 
orientation: 
Disturbed position: 
orientation: 


[-0.15 -0.05 0.39]^, and 

[0.4 0.2 0]^. 

[-0.1499533 -0.050002 0.390070]^, and 

[0.399733 0.199992 3.579856]^. 


The positions of the poles in s-plane are shown in Fig. (5.2). Proportional and derivative 
gain matrices are given below. 
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Figure 5.2: Pole Positions in s-plane for Free Vibration Control of Test Manipulator I 
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619.991 

1311.487 

1598.334 

-101.713 

-1475.722 

885.982 



-22.659 

-883.634 

-1573.340 

-0.690 

1258.533 

-888.109 

Kp = 

= 10® X 

-220.618 

208.040 

561.084 

-47.985 

-443:240 

395.926 


207.997 

595.413 

-139.676 

-21.783 

-389.184 

105.009 



1616.698 

-298.718 

737.603 

-199.752 

-304.517 

-66.533 



_ -2214.663 

-828.419 

-2041.813 

317.959 

1404.277 

-422.330 _ 



862.525 

140.704 

318.976 

-193.860 

-156.365 

677.732 



3.653 

-1154.190 

-2656.574 

80.863 

1264.901 

-710.163 

K„ = 

il 

O 

CO 

X 

-251.774 

-36.448 

1484.511 

18.353 

-469.573 

401.079 


333.842 

1318.403 

-298.774 

4.062 

-348.514 

-68.029 



2556.786 

-755.532 

1861.394 

-400.705 

-160.399 

-112.436 



-4215.574 

-706.665 

-3807.508 

598.105 

1303.164 

-146.496 


The control plots are presented in Fig. (5.3). 


2. Case 2: Test manipulator II 


Desired position; 
orientation: 
Disturbed position: 
orientation; 


[0 0 0.085]^, and 
[0 0 0 ]^. 

(-1.95 X 10"® 6.164 X 10'® 0.0849996]^, and 
[284.762 X 10~® 164.378 x 10"® 2.017 x 10"®]^. 


The positions of the poles in s-plane are shown in the Fig. (5.4). Proportional and 
derivative gain matrices are given below. 


-946.586 

1851.341 

-329.509 

-440.342 

-261.999 

2663.938 

1145.126 ■ 

-1376.650 ■ 

-1029.038 

271.152 

227.342 

-1814.607 

-726.601 ■ 

-1934.165 ■ 

-1097.727 

359.141 

243.555 

-2356.249 

353.132 

2056.321 

-310.414 

-474.549 

-363.133 

2856.158 

1669.7535 

133.309 

-726.146 

24.127 

1.372 

55.613 

-1493.797 

-656.915 

-935.943 

228.617 

128.284 

-1195.752 

-875.909 

1893.424 

197.206. 

-136.271 

-279.742 

1177.966 ' 

147.533 

-1315.242 

-85.325 

93.972 

192.657 

-815.661 

-946.363 

-1835.209 

-33.377 

122.547 

249.190 

-998.750 

228.969 

2177.518 

129.951 

-147.610 

-306.197 

1223.562 

1818.630 

-43.312 

79.275 

-4.164 

5.018 

46.792 

-1699.681 

-849.858 

-58.471 

62.857 

115.339 

-543.643 


The control plots are presented in Fig. (5.5). 
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Translation in Y direction Rotation about Y axis 



Translation in Z direction Rotation about Z axis 


Figure 5.3: Free Vibration Control of Test Manipulator I 











Figure 5.4: Pole Positions in s-plane for Free Vibration Control of Test Manipulator II 


5.3.2 Control of Vibration due to Base Motion 

In this section, simulations of vibration control of the top platform from the base motions 
have been done. For the reason described in the previous section, general gain matrices have 
been used for controlling the top platform motion. Here also, all values are in SI units. 

1. Case 1: Test manipulator I 

Desired position: [—0.15 — 0.05 0.39]^, and 

orientation; [0.4 0.2 0]^. 

Base has six degrees of motion. Motion is assumed to be sinusoidal. In X, V and Z 
direction the frequency is 15 Hz and amplitude is 0.001 m. For Roll, Pitch and Yaw the 
frequency is 10 Hz and amplitude is 0.0001 rad. 

The positions of the poles in s-plane are shown in Fig. (5.6). Proportional and derivative 
gain matrices axe given below. 

840.444 
-794.389 
306.262 
214.406 
-194.647 
-333.970 J 
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Kp = 10® X 


491.093 

1.827 

-143.748 

189.758 

1453.537 

-2392.052 


1472.985 

-1.238.164 

231.715 

847.762 

-252.967 

-1038.071 


1329.266 

-1.108.064 

627.870 

-147.649 

819.228 

-1628.026 


-284.528 

191.276 

-24.983 

-94.791 

-201.008 

466.697 


-1026.719 

792.576 

-175.760 

-448.436 

-10.926 

890.805 
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Translation in Z direction 


Rotation about Z axis 


Figure 5.5: Free Vibration Control of Test Manipulator II 
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Figure 5,6: Pole Positions in s-plane for Controlling Vibrations Due to Base Motions (Test Manip- 
ulator I) 


86.252 

167.020 

290.316 

-31.781 

-119.779 

68.196 

0.365 

-137.375 

-240.525 

18.557 

94.679 

-71.371 

-25.177 

6.353 

140.306 

-0.304 

-26.632 

38.282 

33.384 

134.250 

-40.004 

-6.283 

-42.716 

-2.735 

255.678 

-65.888 

185.377 

-39,345 

4.996 

-15.752 

-421.557 

-98.311 

-358.182 

69.871 

90.135 

-12.013 


The control plots are presented in Fig. (5.7). 
2. Case 2: Test manipulator II 


Desired position; 
orientation: 


[0 0 0.085]^, and 

[0 0 0 ]^. 


The gains are calculated based only on the information that the base platform will be 
subjected to maximum positional disturbance of 0.01 m. No assumption has been made 
regarding nature of the excitation or the magnitude of velocity and acceleration of the 
base due to disturbances. 

The positions of the poles in s-plane are shown in Fig. (5.8). Proportional and derivative 
gain matrices are given below. 
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Figure 5.8: Pole Positions in s-plane for Controlling Vibrations Due to Base Motions of Test Ma- 
nipulator I 


Kp = 10^ X 


-916.842 

1529.756 

-973.280 

246.350 

1886.245 

-1771.092 


1343.907 

-945.913 

-1338.406 

1578.074 

-35.652 

-599.658 


648.639 

-416.583 

-489.040 

653.746 

64.239 

-272.639 


-1548.207 

1058.049 

1385.735 

-1703.212 

-37.787 

682.838 


-0.868.133 

596.751 

781.825 

-962.267 

-3.754 

368.920 


3356.666 

-2305.996 

-2.928.540 

3.612.117 

112.993 

-1494.026 


K„ = 10^ X 


-96.735 172.142 30.323 

161.927 -121.548 -15.324 

-103.035 -169.928 -10.453 
25.870 200.549 23.003 

199.361 -4.328 9.463 

-187.268 -77.528 -10.251 


-61.375 

-52.014 

147.530 

41.975 

35.864 

-102.282 

55.345 

46.110 

-124.011 

-67.823 

-56.944 

153.461 

-1.449 

-0.108 

7.262 

27.052 

21.761 

-67.074 


The base platform is subjected to random excitation in ^-direction as shown in Fig. (5.9). 


The controlled response of the top platform is presented in Fig. (5.10). Some portion of 
Fig. (5.10) is zoomed in Fig. (5.11) to show the response of the top platform clearly. 


5.4 Conclusion 

The simulations shows that even though no particular nature of the disturbances axe assumed, 
the control strategy, developed here, worked quite well. Gains, calculated from the linearized 
dynamic model, should be increased a little at the time of simulating or controlling the actual 
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Figure 5.9: Base Excitation 

system. This control strategy can be used for vibration isolation in many applications e.g. 
isolating telescopes from surrounding disturbances, maintaining the desired direction of the 
naval missiles, etc. After the simulations give satisfactory results, then, the controller should 
be fabricated to control the distmbances. Since the flexibilities of legs are considered here, 
the numerical simulations are closer to the physical system. 
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Chapter 6 
Closure 

6.1 Summary 

In this thesis, singularity measures of a manipulator, dynamics and vibrations of the Stewart 
platform have been investigated. Finally, a vibration control strategy has been developed 
using 6-UPS Stewart platforms. Numerical simulations have been performed to evaluate the 
controller’s performance. 

In chapter 1, a comparison between parallel and serial manipulators and a brief description 
of Stewart platform have been presented. Then the state of the arts have been discussed. 

In chapter 2, a dynamic stability index is introduced and compared with conventional 
measure of static singularity. Though this is illustrated with 6-UPS Stewart platform, it can 
also be used for other parallel manipulators. 

In chapter 3, considering the motion of the base platform and flexibility of the legs, a 
generalized closed form dynamic equation of 6-UPS Stewart platform has been developed 
following the Newton-Euler approach. Using this equations a vibration study of Stewart 
platform has been performed in chapter 4. 

In chapter 5, an active vibration control strategy using Stewart platform has been pre- 
sented. A linearized model has also been developed to calculate the gains of PD controller. 
Finally, simulations of vibration isolation have been done using data of some standard Stewart 
platforms. 
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6.2 Suggestions for Future Work 

The work of this thesis can be extended in many directions. A few among them which may 
have quite significant outcomes, are discussed below. 

1. In chapter 2 significances of the lowest natural frequency of manipulators have been 
discussed. This can be used at the time of designing a manipulator so that manipulators, 
used for a particular pose have a high value of the lowest natural frequency at the desired 
configuration and manipulators designed for tracking purpose should have a high value 
of the lowest natural frequency over the whole path in an overall sense. 

2. The dynamic equation (Eqn. 3.48), derived in the chapter 3, can be non-dimensionalized. 
It will reduce some parameters and also the behaviour of the platform can be predicted in 
terms of non-dimensional parameters. Hence, comparison between two different designs 
of Stewart platform will be much more convenient and accurate. Moreover, one analysis 
will apply to a class of Stewart platform. Scale up or scale down version of a model will 
not need separate analysis. 

3. The generality of the formulation of chapter 3 can be further extended by incorporating 
distribution of the stiffness over the whole leg. Use of Finite Element and Finite Differ- 
ence method will be quite useful in this case. The computation overhead of the FE and 
FD models can be handled by using parallel programming 

4. In chapter 4, vibration study has been done up to httle above the 2nd natural frequency. 
As the simulations have been done in MATLAB, requirements of computation time is 
quite high. Using lower level programming languages (C, C++, FORTRAN, etc.) with 
parallel programming, this impediment can be removed and the work can be extended 
to find out response of the top platform at the base excitation frequency above the 
second natural frequency. Thus, Fig. (4.3) can be completed. The behaviour of the top 
platform at the natural frequencies also can be studied in more details. 

5. In chapter 5, vibration isolation has been done using PD controller. For trajectory 
tracking purpose feedforward strategy can be used. Besides, other control strategies i.e. 
PID, Neural Network, Fuzzy logic etc., can be explored in the future work in this field. 
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Appendix A 

Description of the Test Manipulators 


A.l Test Manipulator I 


All data are in SI unit. 
Base points: 


[ bi b2 ba b4 bs be ] 
Platform points (in platform frame): 

[ Pi P2 P3 P4 Ps Pe ] 


0.6 

0.1 

-0.3 

-0.3 

0.20 

0.5 

0.2 

0.5 

0.3 

-0.4 

-0.30 

-0.2 

0.0 

0.1 

0.1 

0.0 

-0.05 

0.0 


0.3 

0.3 

0.0 

-0.2 

-0.15 

0.15 

0.0 

0.2 

0.3 

0.1 

-0.20 

-0.15 

0.1 

0.0 

0.0 

-0.1 

-0.05 

-0.05 


Unit vectors along fixed axes of imiversal joints: 


[ ki ka ka k4 ks ke ] 


-0.8141 0.2308 
0.2714 0.9231 

0.0000 0.3077 


0.9535 1.0000 0.7071 
0.2860 0.0000 0.7071 
0.0953 0.0000 0.0000 


Mass of lower and upper part of each leg: 


rud = 3.0 and rriu = 1.0 


Centers of gravity of lower and upper parts of each leg (in local frames): 

^dc = [ 0-4 0.14 -0.18 and r^^^ = [ -6.0 —0.08 0.08 ] 


Moments of inertia of lower and upper parts of each leg (in local frames): 



■ 0.010 0.005 0.007 ■ 


I<io = 

0.005 0.002 0.003 
0.007 0.003 0.001 

and I,io = 


0.005 

0.002 

0.002 


0.002 

0.002 

0.001 


0.002 

0.001 

0.003 


-0.9535 

0.2860 

-0.0953 


53 



Centre of gravity of the platform and payload (in platform frame): 

Ro=[ 0.04 0.03 -0.06 ]^ 

Moment of inertia of platform and payload (in platform frame): 

■ 0.050' 0.003 0.004 ' 

Ip= 0.003 0.040 0.003 

[ 0.004 0.003 0.100 _ 

Coefficients of viscous friction: 

Cu = 0 . 0001 , Cp = 0 . 001 , Cs = 0.0002 

Stiffness of legs: 


Klegi 

= 2 

X 

10® 

^leg2 

= 2 

X 

10® 

^legz 

= 2 

X 

10® 

^leg^ 

= 2 

X 

10® 

^legs 

= 2 

X 

10® 


= 2 

X 

1 n8 

XU 


A. 2 Test Manipulator II 

All data are in SI unit. 

Base points: 

■ 0.03 0.004486 -0.015 -0.02793 -0.015 0.02345 

[ bi b2 bs b4 bs be ] = 0.00 0.02966 0.02598 -0.01095 -0.02598 -0.01872 

[ 0.00 0.0 0.0 0.00 -0.0 0.0 

Platform points (in platform frame): 

■ 0.0428 0.0326 -0.0334 -0.0432 -0.0094 0.0106 

[ Pi P2 P3 P4 P5 P6 ] = 0.0139 0.0313 0.0301 0.0127 -0.0440 -0.0437 

• [ 0.0 0.0 0.0 0.0 0.0 0.0 
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Mass of lower and upper part of each leg: 

rUd = 0.3 and rriu = 0.1 

Centers of gravity of lower and upper parts of each leg (in local frames): 

Tdo = [ 0.04 0.014 -0.018 f and = [ -0.6 -0.008 0.008 

Moments of inertia of lower and upper parts of each leg (in local frames): 

■ 0.000010 0.000005 0.000007 1 [ 0.000005 0.000002 0.000002 ‘ 

Ido = 0.000005 0.000002 0.000003 and = 0.000002 0.000002 0.000001 

[ 0.000007 0.000003 0.000001 J [ 0.000002 0.001 0.003 

Centre of gravity of the platform and payload (in platform frame); 

Ro = [ 0 0 0 

Moment of inertia of platform and payload (in platform frame); 

■ 0.0001 0.0 0.0 

Ip = 0.0 0.0001 0.0 

[ 0.0 0.0 0.0001 _ 

Coefficients of viscous friction: 

Cu = 0.0001, Cp = 0.001, a = 0.0002 

Stiffness of legs: 


^legi 

= 8 X 10^ 

Kleg2 

= 8 X 10^ 

Kleg2 

= 8 X 10^ 

^ItgA 

= 8 X 10'^ 

^legs 

= 8 X 10^ 

^lege 

= 8 X 10^ 


All other parameters are same as of previous section. 
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Appendix B 

Jacobian Expression 


Jacobian expression (A), required for devloping the linearized model of the chapter 5, is 
developed here. 

A(1:6,1:6) = 0 (B.l) 

A(1 : 6,7 : 12) = [6 x 6] Identity Matrix (B.2) 

Let, 

a = J-(HF+[^f“]-, + d) (B.3) 

and. 
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where 


^ ^ d{HF) d\ 3?pFext 

dzi dzi dzi dzi _ 3?pMe3.4 

aJp, ^ ajpiat 

dzi dzi dzi ^ 

i=i 


Jpiat = Moment of inertia of top platform. 
Jpj = Moment of inertia of :?’-th leg. 


(B.8) 

(B.9) 


dzi 


^Jplat 

dzi 

dME^ 

dzi 


dMEs 

dz^ 

aMRp 


dMKp 

dzi 


^{I + M(i?2E3-RpR^)} J 


= 0 


aMR 

dzi 


= M^SRpRo 


0 for z 7 ^ 4, 5 or 6 


(B.IO) 

(B.ll) 

(B.12) 


dz4 

a3?p 

dzs 

dzQ 


0 c{zs)s(z5)c{z^) + s{z6)s{z4) - c{ze)s{z5)s{z4) + s{z6)c{z^) ' 

0 siz6)s{z5)c{z4) - c(z6)s(z4) -s(ze)s(z5)s(z4) - c{z6)c{z4) 

0 c(z5)c{z4) -c{z5)s{z4) 

-c(zg)s{z5) c(z6)c{z5)s{z4) c{ze)c{zs)c{zA 

-s(z6)s{z5) s(z6)c{z5)s{z4) 3(ze)c(zs)c(z4) 

-c(zs) -s(zs)s(z4) -s(zs)c(z4) 

-s(ze)c(zg) -s(Ze)s(Z5)s(Z4) - c(ze)c(z4) -s(z6)s(z5)c(z4) + c(z6)s(z4) 
c(ze)c(zs) c(z6)s(z5)s(z4) - s(ze)c(z4) c(ze)s(z5)c(z4) + 5 ( 26 ) 5 ( 2 ^) 

0 0 0 


Where, c(.) = cos(.) and s(.) = sin(.). 


The (2,2) element of Jput [in Eqn. (B.IO)] is 


^{I + M(i?jE3-RpRj)} 

OZi 




dzi 


dzi 


(B-13) 


-.dUr. 


a3R„ 


-R,R5.(^n. 


33 

Using expressions (B.ll), (B.12) and (B.13) we can find out the expression of from 

Eqn. (B.IO). 
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Expanding the 2nd term of expression (B.9) we get, 


dzi 


9Qpj g(Qpjqp,) 

dzi dzi ^ 

^(qpjQpj) ^(qp,Qp,qp,) 

dzi dzi 

dQpj 

dzi 


-(?.9a q~ J_Q fSEi') 
V dzi qpj ^ ^Pi dzi > 


- 4 -a~ ^ -r^O 4 -a Q ^ 

dzi ^Pi ^ qpj dzi \ dzi ^pj^pj +qpj dzi qpj + qpj'->ipj g^. . 

Other than and all terms of the Eqn. (B.14) axe known. 

dqj, . dRp 


(B.14) 


dzi dzi 


At the time of developing expressions for the suffixes (i and j) have been dropped to 
reduce the cluttering of the expressions. Before going into derivation of those expressions 
which axe required for derivation have been evaluated below. 


Pi 


(B.15) 


dz 

ds 

Yz 

all k X s II 






dz dz 


d_L, 

dz' 


1 , as ^ 


as, 


dz 


2 I k X s 


{(kx^)^(kxs) + (kxs)^(kxg)} 


dz 


dT 

dz 


_a 

dz 


k X s s x (k x s) 


k X s 


k X s 


^{s(Irf + I„)s} = |i(I, 4 -I„)s+s(|^ + ^)s + s(Id + Ij|j 

We can break in the four following terms. 

aQp 


dz 

Where 

Terml 


- = Terml + Term2 — Term3 — Term4 


(B.17) 

(B.18) 

(B.19) 

(B.20) 

(B.21) 


[2m, 


(t.r„ + s.||)T- f (s.r,) 


+ mu((^)^ru + r^^)}L^ - 2L^{mdr\ + 
iss 


+K(i + + s(|)n 

Term 2 = ijE 3 lW(^rr. + rj^) + 


(B.22) 
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Terms 


(B.23) 

(B.24) 


-2T^(mrfr^ + rriurl)] 


dL 

dz 


(sr^ + r„s^)] 


Term4 = ^[L^{md[{^^ x r^ + s x ^)(s x Vdf + (s x rd)(^ x + s x ^7] 

:ds , ,,ds 

V r 4- fi y 

^dz 

4-^[s(Id + Iu)s]} - {md{s X r<i)(s x + m„(s x r„)(s x 


r ^ KytJ KJ M. >tl ^ , ^rp . . . lyo dru.'T 

+’^u[(-^ X r„ + s X -^)(s X r„)^ + (s x r„)(— x r„ + s x 


dz 


)"] 


+s(Id + I„)s } 2 X 


dL 

dz^ 


(B.25) 


^Qp . 


Using Eqn. (B.16), (B.17), (B.19), (B.20) and Eqn. (B.22) to (B.25) we can express 




in terms of platform parameters and z. Again using Eqn. (B.15) we can evaluate Hence, 


dJ 

this completes the expression of in' terms of platform parameters and z. 

Next we are going to find out the expression of For that we require the expressions of 


Qhf ^ _a. 

dzi ’ dzi 8z. 


dzi 


^pFext 

3^h4ea;t 


where 


SHF _ dH dF 
dzi 'dzi'^^dzi 


(B.26) 


& 

V ''j-th column 


dzi 


dzi 

/ ^Qpj . ^ dsj \ 

(■^ ^ ®i + %■ ^ dzi) 


K ^ 

^^9} dzi 


dziJ J 


jf = 1 to 6 ; 


(B.27) 


(B.28) 


1x6 


dr] 

dzi 

drjpiat 

dZi 


9Vj 

dzi 


^Vplat . ^ 


dzi dz. 


) ” 

rE% 


i=i 


X (a,, X R,) + w, X X R,) + Wp X (Wp X ^R<,,)} 

^xH + u;,xga,p + arpXl^ 

+M(^R<i,)x{(«,-RK-g} 

+MRp X {(^ ■ RpVp + {wp - (^R<ipVp) + (up ■ Rp)^} 

dVi 


dzi 

^xV +a , 

dzi ^ ^ dZi dZi J 


(B.29) 


(B.30) 


(B.31) 
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Expression of has been developed later. Other terms of Eqn. (B.29), (B.30) and (B.31) 
have already been found out. 


dzi 


3^Fext 


. 1 

dzi ^ 

_ 3?pMext . 


. ^Mext _ 


(B.32) 


Expression of ^ can be found out using (B.26), (B.29) and (B.32). 

Hence we have got the expression for Jacobian matrix (A), in terms of Stewart platform 
parameters and z. 


as. 

dzi 

dWj 

dzi 

dij 

dzi 

aui, 

dzi 

dzi 


aus, 

dZi 


aus, 

dzi 


dzi 


aus, 

dzi 


dojp dcu atp. 

^{L,(|xS,+s,xf)-||(s,xS.)} 

asj 


duj. 


dSj 

dzi 


/ \ , /9u}p , . aqp. . 

^ X (wp X qpj d-Wp X (-^ X qp,) +a;p X {Up X -^) 

+(S, - L,s,n^ - |is, - - (S, - 4s,)>|il 

1 aui, ^,aL,„, -aw,-., 

g|{ g;: X U, + s, X ^ - 2(^w, + L,^)}L, 


a^i 

X U,, - 24W,.)1 


au 


^xr,,+U.,x.^ + -g^x(W,xr4 

dUj dSj 8V2, ,, ar„, aw, , 

"ST &r “a?" 

■ - -- 


X aw, 


^xr,) + W,x(W,x^) + 2(||w,xs, 

ywj dsj. 

^xs,+L,W,x_) 

arrf. aus ar„ ^^4 

d aif ■" ^ "Sf' ■" ^ 

aid, ai„,,„ „ , ,au2, aw,- ,, , 


(B.33) 

(B.34) 

(B.35) 

(B.36) 

(B.37) 

(B.38) 

(B.39) 


(B.40) 
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C)c. 

{^nji-^.Vi. + S 


1 '•"•'t 

, r* ^s,- . 1 Qg. 

dz, ’■^^'’'^“"*“'S:-s}"-I?{%(^xU5, 


(B.41) 


+s 


dUs. dL- ' 

~ ^ Us,)} + (m„,s,-.U4, + C„-4 - m„,s,.g)|!i (B.42) 

UZi 
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